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I The notions of generating sets of conservation laws of systems of differential equations with 

respect to symmetry groups and equivalence groups are introduced and applied. This 
\ allows us to generalize essentially the procedure of finding potential symmetries for the 

systems with multidimensional spaces of conservation laws. A class of variable coeffi- 
cient (l-l-l)-dimensional nonlinear diffusion-convection equations of general form f{x)ut = 
. {g{x)A{u)ux)x + h{x)B[u)ux is investigated. Using the most direct method, we carry out 

two classifications of local conservation laws up to equivalence relations generated by both 
t~i . usual and enhanced equivalence groups. Equivalence with respect to G~ and correct choice 

I of gauge coefficients of equations play the major role for simple and clear formulation of the 

final results. The notion of contractions of conservation laws and one of characteristics of con- 
servation laws are introduced and contractions of conservation laws of diffusion-convection 
equations are found. 

> 

m 

^ ■ 1 Introduction 
cn 

The presented paper continues the series of works on modern trends of group analysis illustrated 
on example of nonlinear variable coefficient diffusion-convection equations 

o 

O ■ f{x)ut = {g{x)A{u)ux)x + h{x)B{u)ux (1) 



started in [14,15]. Here / = /(x), g = g{x), h = h{x), A = A(u) and B = B{u) are arbitrary 
smooth functions of their variables, f{x)g{x)A{u)y^O. 
5^ I In the first two parts of the presented series we have shown an importance of detailed study 

of different kinds of equivalence transformations for solving the group classification problem and 
some its applications. In particular, we investigated equivalence groups (usual and extended 
ones) , discussed their structure and performed the complete group classification of class ([T]) . We 
determine that the complete usual equivalence group G~ of class ([1]) consists of the transforma- 
tions 



i = 6it + 62, X = X{x), u = 63U + 64, 

f = ^^f, g = eie2^X^g, h = eie^^h, A = e2A, B = e3B, 

where 6j {j = 1,4) and £i {i = 1,3) are arbitrary constants, 5i53eie2£3 7^ 0, X is an arbitrary 
smooth function of x, 7^ 0. 

It appears that class ([1]) admits further equivalence transformations which do not belong 
to G~ and depend on arbitrary elements in some non-fixed (possibly, nonlocal) way. We have 
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constructed the complete in this sense extended equivalence group G"" of class ([T]), using the 
direct method. G"" is formed by the transformations 

i = 6it + 62, x = X{x), u = 53U + 54, 

S S 

f = Arr^f, g = eie^^X^ipg, h = eiel^iph, A = e2A, ^ = £3(5 + 64^), 

where 6j (j = 1,4) and £i {i = 1,4) are arbitrary constants, 5i53eie2£3 7^ 0, X is an arbitrary 

smooth function of x, 7^ 0, 99 = e <^ . 

It is shown that application of extended equivalence group G~ is of vital importance for 
obtaining group classification in closed explicit form [13,14]. Likewise we state below the ne- 
cessity of application of equivalence transformations to investigation of conservation laws. More 
precisely, it will be shown that classification of conservation laws of equations with respect to 
the usual equivalence group C can be formulated in an implicit form only. At the same time, 
using the extended equivalence group G^, we can present the result of classification in a closed 
and simple form with a smaller number of inequivalent equations having nontrivial conservation 
laws. 

Using the gauge equivalence transformation t = t, x = f ^jy, u = u we can reduce equa- 
tion ([H) to 

f{x)u^ = {A{u)us;)i + h{x)B{u)ui, 

where /(x) = g{x)f{x), g{x) = 1 and h{x) = h[x). That is why, without loss of generality we 
can restrict ourselves to investigation of the equation 

f{x)ut = {A{u)ux)^ + h{x)B{u)ux. (2) 

Any transformation from which preserves the condition 5 = 1, has the form 

i=5it + S2, x = 5^j e^^ f^dx + 6e, u = 63U + 6^, 

f = SiS^H^fe-^^^ /^ h = 6g6:f'he-^^ /^ (3) 
A = 655gA, B = 67{B + 5sA), 

where 6i {i = 1,9) are arbitrary constants, Sidsd^djdg 7^ 0. (Here and below jh = Jh{x)dx.) 
The set of such transformations is a subgroup of G^ . 

Result of group classification allows us to construct a number of exact solutions of differ- 
ent equations from class ([1]). In the second part [15] of the series we have used two different 
approaches for finding exact solutions: direct application of Lie (conditional and generalized 
conditional) symmetries and reconstruction of new solutions by acting on the known ones with 
additional equivalence transformations. 

After analyzing the obtained results it appears that some of the previously classified equations 
can be regarded as limiting cases of the other ones, that leads to the natural notion of contractions 
of (systems of) equations introduced in the second part [15] of the series. Ibid we investigated 
contractions of equations from class ([1]). 

These ideas can be generalized for investigation of conservation laws of systems of differential 
equations. Thus, e.g., there exists a close similarity between characteristics of conservation 
laws and symmetries that leads to the natural notions of contractions of characteristics and 
conservation laws that will be considered below. These notions are illustrated by example of 
conservation laws of equations ([T]). 

Similarly to the group classification problem, we show that the complete description of the 
spaces of conservation laws of equations ([I]) is impossible without essential usage of (usual and 
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extended) equivalence transformations. Investigation of such equivalences allows us to simplify 
essentially technical calculations and to obtain results in closed explicit form. Equivalence with 
respect to the extended equivalence group and correct choice of gauge coefficients of equations 
allow us to obtain clear formulation of the final results. Note that for wide classes of nonlinear 
equations, this is the only way to obtain the complete description of spaces of conservation 
laws (see, e.g., [12,13,26] for detailed analysis of application of such equivalences for finding 
conservation laws of different classes of diffusion-type equations). 

The notion of equivalence of conservation laws with respect to a group of transformations 
which was introduced in [26] can be generalized in several directions: classification of pairs 
"system + space of conservation laws"; classification of conservation laws for a given system 
with respect to its symmetry group; classification of pairs "system + a conservation law". Such 
classification schemes are helpful both for ordering and deeper understanding of the found con- 
servation laws and for further applications. For instance, investigation of different generating 
sets of conservation laws of equations ^ allows us to close the "blank spot" in theory of poten- 
tial symmetries. Namely, previously, for construction of simplest potential systems in cases when 
the dimension of the space of local conservation laws is greater then one, only basis conservation 
laws were used. However, the basis conservation laws may be equivalent with respect to groups 
of symmetry transformations, or vice versa, the number of C-independent linear combinations 
of conservation laws may be grater then dimension of the space of conservation laws. The first 
possibility leads to an unnecessary, often cumbersome, investigation of equivalent systems, the 
second one makes possible missing a great number of inequivalent potential systems. Below 
we show how to choose conservation laws in order to obtain all possible inequivalent potential 
systems associated to the given system. 

In Section [2] we adduce basic definitions and statements on conservation laws. An impor- 
tant notion of characteristics of conservation laws is discussed in Section [3l Following the spirit 
of [26] we repeat a definition of equivalence of conservation laws with respect to a transformation 
group in Section HI In Section S] we introduce also the notions of generating sets of conservation 
laws with respect to symmetry groups and with respect to equivalence groups. After that (Sec- 
tion [5|) two different classifications of conservation laws of equations ([1]) are presented, namely, 
classification with respect to the usual equivalence group and one with respect to the extended 
equivalence group. Next (Section [6]), we introduce the notions of contractions conservation laws 
and ones of characteristics of conservation laws and give some examples of the contractions in 
class (dJ . Generating sets of conservation laws of diffusion-convection equations ([1]) are investi- 
gated in Section [71 Such investigation form a basis for finding all possible inequivalent potential 
systems of equations ([1]) presented in Section[9l In Section[8]we discuss shortly potential symme- 
tries and potential conservation laws in general case and in case of systems with two independent 
variables. Potential conservation laws of equations ([1]) are investigated in Section [TOl 

A detailed investigation of potential symmetries obtained from the constructed potential 
systems will be performed in the last part [16] of this series. 

2 Basic definitions and statements on conservation laws 

In this and next two sections we give basic definitions and statements on conservation laws and 
formulate the notion of equivalence of conservation laws with respect to equivalence groups, 
which was first introduced in [26] and some of its generalizations. This notion is a basis for 
modification of the direct method of construction of conservation laws, which is applied in 
Section [5] for exhaustive classification of local conservation laws equations from class ([1]) . 

Let £ be a system L{x, U(p)) = of / differential equations = 0, . . . , L' = for m unknown 
functions u = (n^, . . . , u"*) of n independent variables x = {xi, . . . , x„). Here ii(p) denotes the set 
of all the derivatives of the functions u with respect to x of order no greater than p, including 
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u as the derivatives of the zero order. Let /^(^^ denote the set of all algebraically independent 
differential consequences that have, as differential equations, orders not greater than k. We 
identify with the manifold determined by in the jet space 

Definition 1. A conserved vector of the system C is an n-tuple F = {F^{x, • • • , F"'{x, 
for which the divergence DivF = Dj-F* vanishes for all solutions of C (i.e. DivF|^= 0). 

In Definition [1] and below Di = denotes the operator of total differentiation with respect 
to the variable Xj, i.e., Di = 5^. + u'^idu^, where u'^ and u'^^ stand for the variables in jet 

spaces, which correspond to derivatives d^°'^u"' /dx"^ . . . dx"" and du'^/dxi, a = {ai, . . . ,an), 
€ N U {0}, \a\: = ai + • • • + a^. We use the summation convention for repeated indices and 
assume any function as its zero-order derivative. The notation means that values of V are 
considered only on solutions of the system C. 

Definition 2. A conserved vector F is called trivial if F'^ = F^ A- F', i = l,n, where and 
F* are, likewise F*, functions of x and derivatives of u (i.e., differential functions), F"* vanish 
on the solutions of C and the n-tuple F = (F^, . . . , F") is a null divergence (i.e., its divergence 
vanishes identically). 

The triviality concerning the vanishing conserved vectors on solutions of the system can 
be easily eliminated by confining on the manifold of the system, taking into account all its 
necessary differential consequences. A characterization of all null divergences is given by the 
following lemma (see e.g. [23]). 

Lemma 1. The n-tuple F = (F^, . . . ,F"), n>2, is a null divergence (DivF = 0) iff there exist 
smooth functions v^^ (i,j = l,n) of x and derivatives of u, such that v^^ = —v^^ and F* = Djv'^K 

The functions w*-' are called potentials corresponding to the null divergence F. If n = 1 any 
null divergence is constant. 

Definition 3. Two conserved vectors F and F' are called equivalent if the vector- function F' — F 
is a trivial conserved vector. 

The above definitions of triviality and equivalence of conserved vectors are natural in view 
of the usual "empiric" definition of conservation laws of a system of differential equations as 
divergences of its conserved vectors, i.e., divergence expressions which vanish for all solutions of 
this system. For example, equivalent conserved vectors correspond to the same conservation law. 
It allows us to formulate the definition of conservation law in a rigorous style (see, e.g., [28]). 
Namely, for any system C of differential equations the set CV{C) of its conserved vectors is a 
linear space, and the subset CVo(£) of trivial conserved vectors is a linear subspace in CV{C). 
The factor space Ch{C) = CV(/3)/ C\q{C) coincides with the set of equivalence classes of CV(£) 
with respect to the equivalence relation adduced in Definition [3l 

Definition 4. The elements of CL(jC) are called conservation laws of the system C, and the 
whole factor space CL(/^) is called the space of conservation laws of C. 

That is why we assume description of the set of conservation laws as finding CL(£) which 
is equivalent to construction of either a basis if dimCL(>C) < oo or a system of generatrices 
in the infinite dimensional case. The elements of CV(£) which belong to the same equivalence 
class giving a conservation law are considered all as conserved vectors of this conservation 
law, and we will additionally identify elements from CL(i2) with their representatives in CV(£). 
For F E CV(£) and J- G CL(£) the notation F £ will denote that F" is a conserved vector 
corresponding to the conservation law J^. In contrast to the order rp of a conserved vector F as 
the maximal order of derivatives explicitly appearing in F, the order of the conservation law T 
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is called mm{rp \ F S J^}. Under linear dependence of conservation laws we understand linear 
dependence of them as elements of CL(£). Therefore, in the framework of "representative" 
approach conservation laws of a system C are considered linearly dependent if there exists a 
linear combination of their representatives, which is a trivial conserved vector. 

3 Characteristics of conservation laws 

Let the system C be totally nondegenerate [23] . Then application of the Hadamard lemma to the 
definition of conserved vector and integrating by parts imply that divergence of any conserved 
vector of C can be always presented, up to the equivalence relation of conserved vectors, as a 
linear combination of left hand sides of independent equations from £. with coefficients A'^ being 
functions on a suitable jet space 

j(k). 



Here the order k is determined by C and the allowable order of conservation laws, /i = 1,1. 

Definition 5. Formula ^ and the /-tuple A = (A^, . . . , A') are called the characteristic form 
and the characteristic of the conservation law Div F = correspondingly. 

The characteristic A is trivial if it vanishes for all solutions of C. Since C is nondegenerate, 
the characteristics A and A satisfy ^ for the same F and, therefore, are called equivalent iff 
A — A is a trivial characteristic. Similarly to conserved vectors, the set Ch(£) of characteristics 
corresponding to conservation laws of the system £ is a linear space, and the subset Cho(£) of 
trivial characteristics is a linear subspace in Ch(£). The factor space Chf(£) = Ch(£)/ Cho(£) 
coincides with the set of equivalence classes of Ch(£) with respect to the above characteristic 
equivalence relation. 

The following result [23] forms the cornerstone for the methods of studying conservation laws, 
which are based on formula ([4]), including the Noether theorem and the direct method in the 
version by Anco and Bluman [1,2]. 

Tlieorem 1 ( [23]). Let C be a normal, totally nondegenerate system of differential equations. 
Then representation of conservation laws of C in the characteristic form ^ generates a one- 
to-one linear mapping between CL(£) and Chf(£). 

Using properties of total divergences, we can exclude the conserved vector F from ^ and 
obtain a condition for the characteristic A only. Namely, a differential function / is a total 
divergence, i.e., / = DivF for some n-tuple F of differential functions iff E(/) = 0. Hereafter 
the Euler operator E = (E^, • • • , E™) is the m-tuple of differential operators 



E« = (-Z))-5„., a = l,m, 

where a = (ai, . . . , ««) runs the multi-indices set (oj gNU{0}), (—1))" = (— Di)"^ . . . {—Dm) 
Therefore, action of the Euler operator on Q results to the equation 



which is a necessary and sufficient condition on characteristics of conservation laws for the sys- 
tem C The matrix differential operators and D2 ^-re the adjoints of the Frechet derivatives D;^ 
and D^,, i-e., 



DivF = A'^L'^. 



(4) 



(5) 
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Since D^(-^^) = automatically on solutions of L then equation ([5]) implies a necessary condition 
for A to belong to Ch(£): 

Di(A)|^=0. (6) 

Condition ([6]) can be considered as adjoint to the criterion D/^(?7)|£= for infinitesimal invariance 
of L with respect to evolutionary vector field having the characteristic r\ = {rj^, . . . ,7]"^)- That 
is why solutions of ([6]) are called sometimes as cosymmetries [4] or adjoint symmetries [2]. 



4 Equivalence of conservation laws 

with respect to transformation groups 

We can essentially simplify and order classification of conservation laws, taking into account 
additionally symmetry transformations of a system or equivalence transformations of a whole 
class of systems. Such problem is similar to one of group classification of differential equations. 

Proposition 1. Any point transformation g maps a class of equations in the conserved form 
into itself. More exactly, the transformation g: x = Xg{x,u), u = Ug{x,u) prolonged to the 
jet space J^^'^ transforms the equation DiF^ = to the equation DiF^ = 0. The transformed 
conserved vector Fg is determined by the formula 

D . X ■ 1 

Fgix,U(r)) = i^'-j F^{x,U^r)), ^-C Fg{x,U^r)) = JjjY\iDxX)F{x,U(^r)) (7) 

in the matrix notions. Here \Dx x\ is the determinant of the matrix F)xX — (Z^^ x^). 

Note 1. In the case of one dependent variable (m = \) g can be a contact transformation: 
X = Xg(x,U(i)), -£((1) = iig(i)(x, U(i)). Similar notes are also true for the statements below. 

Definition 6. Let G be a symmetry group of the system C. Two conservation laws with the 
conserved vectors F and F' are called G- equivalent if there exists a transformation g £ G such 
that the conserved vectors Fg and F' are equivalent in the sense of Definition O 

Any transformation g £ G induces a linear one-to-one mapping g^ in CV(>C), transforms 
trivial conserved vectors only to trivial ones (i.e., C\o{C) is invariant with respect to g*) and 
therefore, induces a linear one-to-one mapping g^ in CL(>C). It is obvious that gf preserves linear 
(in)dependence of elements in CL{C) and maps a basis (a set of generatrices) of CL{C) in a 
basis (a set of generatrices) of the same space. In such way we can consider the G-equivalence 
relation of conservation laws as well-determined on CL(>C) and use it to classify conservation 
laws. 

Proposition 2. // the system C admits a one-parameter group of transformations then the 
infinitesimal generator X = + rj^du'^ of this group can be used for construction of new 
conservation laws from known ones. Namely, differentiating equation ^ with respect to the 
parameter e and taking the value e = 0, we obtain the new conserved vector 

r = -X^r)F' + (DjOF^ - {Dji^)F\ (8) 

Here X^^.) denotes the r-th prolongation [23, 24] of the operator X . 

Note 2. Formula ([8]) can be directly extended to generalized symmetry operators (see, for 
example, [5,18]). A similar statement for generalized symmetry operators in evolutionary form 
iC = 0) was known earlier [9,23]. It was used in [19] to introduce a notion of basis of conservation 
laws as a set which generates a whole set of conservation laws with action of generalized symmetry 
operators and operation of linear combination. 
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Proposition 3. Any point transformation g between systems C and C induces a one-to-one 
linear mapping g^, from CV(£) into CV(£), which maps C\q{C) into C\q(C) and generates a 
one-to-one linear mapping g^ from CL(£) to CL(£). 

Corollary 1. Any point transformation g between systems C and C induces a one-to-one linear 
mapping gf from Ch[{C) to Chf(£). 

It is possible to obtain an explicit formula for correspondence between characteristics of C 
and C. Let = Af'^C, where A'"' = A'""^D°', A^""^ are differential functions, a = (ai, . . . , a„) 
runs the multi-indices set (oj €N U {0}), ^, = 1, /. Then 

Here A'^'^* = (— Z))" • A^'^°' is the adjoint to the operator A'^^. For a number of cases, e.g. if C 
and C are single partial differential equations (/ = !), the operators A^'^ are simply differential 
functions (i.e. A^"'" = for |q| > 0) and, therefore, A''^^* = A^" . 

Consider the class C\s of systems Ce- L{x,U(^p^,9{x,U(^p-^)) = parameterized with the 
parameter-functions 9 = 9{x, U(^p-^). Here L is a tuple of fixed functions of x, U(p) and 9. 9 denotes 
the tuple of arbitrary (parametric) functions 9{x, u^p-^) = {9^{x, Uf^p-^), . . . , 9^{x, ui^p-^)) running the 
set S of solutions of the system S{x,U(^p-j,9(g){x,U(p))) = 0. This system consists of differential 
equations on 9, where x and U(p) play the role of independent variables and 9^^^ stands for the 
set of all the partial derivatives of 9 of order no greater than q. In what follows we call the func- 
tions 9 arbitrary elements. Sometimes the set S is additionally constrained by the non- vanish 
condition S' {x,U(^p-^,9(^q^{x,U(^p^)) ^ with another tuple S' of differential functions. Denote the 
point transformations group preserving the form of the systems from C\s as G~ = G'^{L, S). 

Consider the set P = P{L,S) of all pairs each of which consists of a system Cq from C\s 
and a conservation law of this system. In view of Proposition [3l action of symmetry transfor- 
mations of the system Cg and transformations from C on C\s and {CY{Cg) \ 9 £ S} together 
with the pure equivalence relation of conserved vectors naturally generates equivalence relations 
on P. Such generations can be made in several directions: classification of pairs "system -|- 
space of conservation laws" ; classification of conservation laws for a given system with respect 
to its symmetry group; classification of pairs "system -|- a conservation law". Let us consider 
these possibilities in more detail. 

1. Equivalence with respect to an equivalence group [26]. We wish to find a list of 
systems from the class C\s where (i) Cg are graded by the dimensions of conservation laws and 
(ii) we can choose a representative, all other systems are equivalent to, or we can write conditions 
that systems belong to a given class. 

Definition 7. Let 9,9' £ S, T e CL{Ce), J" G CL(/:e'), F ^T,F' £ T' . The pairs {Le^T) 
and {Cgf^J-') are called -equivalent if there exists a transformation g € C which transform 
the system Ce to the system Cqi and such that the conserved vectors Fg and F' are equivalent 
in the sense of Definition [3l 

Classification of conservation laws with respect to G~ will be understood as classification 
in P with respect to the above equivalence relation. This problem can be investigated in a way 
that is similar to group classification in classes of systems of differential equations, especially if 
it is formulated in terms of characteristics. Namely, we construct firstly the conservation laws 
that are defined for all values of the arbitrary elements. (The corresponding conserved vectors 
may depend on the arbitrary elements.) Then we classify, with respect to the equivalence group, 
arbitrary elements for each of that the system admits additional conservation laws. 

In an analogues way we also can introduce equivalence relations on P, which are generated 
by either generalizations of usual equivalence groups or all admissible point or contact transfor- 
mations (called also form-preserving ones [21]) in pairs of equations from C\s- 
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2. Generating sets of conservation laws. Consider system Lq and a non-empty subset 
5^ = {F^, 7 € r} of the set QN(^Lq) of its conserved vectors. Acting on by the transformations 
from the symmetry group G^^^i^Lo) of Lq and taking hnear combinations, we may obtain wider 
set of conserved vectors. Roughly speaking, the set ^ is independent with respect to Q^'^^ if an 
arbitrary nonzero Hnear combination of G^'^^-transformed conservation laws is nontrivial. The 
rigorous definition is as follows. 

Definition 8. The subset 5^ = {Fy,7 G T} C QN(^IIq) is dependent with respect to the 
group G^^^{Cq) of symmetry transformations if 

V{7i, . . . , 7n} C r, G KolF^ = ^ A, = 0, i = T;^^. 

The set of G™^^(>C0)-independent conserved vectors generating the whole space C\{Ce) is 
called generating set for the given system. Since the operations of symmetry transformation and 
taking the linear combination commute, the above definition is well-posed. 

Generating sets of such type, called also L-bases, were firstly described in [19] (see also [9]). 

3. Generating sets of pairs "system -|- conservation law" . Consider now the action of an 
equivalence transformation on a pair from P. If the first element of the resulting pair coincides 
with the first element of another pair from P, we can take linear combinations of their second 
elements. All conservation laws that cannot be obtained in such way will be called inequivalent. 
In such way we obtain a generating set of pairs "system + conservation law" . 

It becomes especially interesting if such generating set depends on arbitrary elements (pa- 
rameters) of the system. In particular, it is possible that the characteristics of conservation laws 
do not depend on the parameters. Then, the action of equivalence group on a pair from the 
generating set maps it to a pair from P. Taking the linear combinations of conservation laws 
for the pairs with coinciding first elements, we get all possible conservation laws. 

Note 3. A notion of generating sets of pairs "system + conservation law" with respect to all 
admissible transformations is introduced in an analogous way. 

Note that a problem of constructing such generating sets is different from the one of classifica- 
tion of conservation laws with respect to equivalence group: under the problem of classification 
of conservation laws with respect to equivalence group we understand classification of pairs 
"equation+ conservation law", and a conservation law from the generating set can generate a 
subspace of the space of conservation laws for equation with another values of arbitrary elements. 

Note 4. It can be easy shown that all the above equivalences are indeed equivalence relations, 
i.e., they have the usual reflexive, symmetric and transitive properties. 

5 Local conservation laws of diffusion— convection equations 

We search (local) conservation laws of equations from class ([2]) , applying the modification of the 
direct method, which was proposed in [26] and is based on using the notion of equivalence of 
conservation laws with respect to a transformation group and classification up to the equivalence 
group of a class of differential equations. 

Conservation laws were investigated for some subclasses of class ([I]). In particular, conser- 
vation laws of the linear heat equation and Burgers equation were constructed firstly in [3]. 
V.A. Dorodnitsyn and S.R. Svirshchevskii [8] (see also [10, Chapter 10]) completely investigated 
the local conservation laws for one-dimensional reaction-diffusion equations ut = {A{u))xx+C{u) 
having non-empty intersection with class ([T]). In [18] the first-order local conservation laws of 
equations ([T]) are found. Developing results obtained in [7] for the case hB = 0, / = 1, in 
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the recent papers [11,26] we completely classified potential conservation laws (including arbi- 
trary order local ones) of equations ([T|) with f = g = h = 1 with respect to the corresponding 
equivalence group. Conservation laws of equations ([T|) were considered in [13]. 

In view of results of Section H] it is sufficient for exhaustive investigation if we classify conser- 
vation laws of equations only from class ([2]). 

There are two independent variables t and x in equations under consideration, which play a 
role of the time and space variables correspondingly. Therefore, the general form of constructed 
conservation laws will be 

DtF{t, X, + DxG{t, X, ■U(r)) = 0, (9) 

where Dt and Dx are the operators of total differentiation with respect to t and x. The com- 
ponents F and G of the conserved vector {F, G) are called the density and the flux of the 
conservation law. Two conserved vectors {F,G) and {F',G') are equivalent if there exist such 
functions F, G and H olt, x and derivatives of u that F and G vanish for all solutions of C and 
F' = F + F + DxH, G' = G + G- Df H. 

At first we prove the lemma on order of local conservation laws for more general class of 
second-order evolution equations, which covers class ([I]). 

Lemma 2. Any local conservation law of any second-order (1 + 1)- dimensional quasi-linear 
evolutionary equation has the first order and, moreover, there exists its conserved vector with 
the density depending at most on t, x, and u and the flux depending at most on t, x, u and Ux- 

Proof. Consider a conservation law of form ([9]) of the second-order (1 + l)-dimensional quasi- 
linear evolutionary equation 

Ut = S{t,X,U,Ux)Uxx + R{t,X,U,Ux), (10) 

where S ^ 0. In view of equation (jlOp and its differential consequences, we can assume that F 
and G depend only on t, x and Uk = d^u/dx^, k = 0,r', where r' < 2r. Suppose that r' > 1. 
We expand the total derivatives in ([9]) and take into account differential consequences of form 
utj = Dx{Suxx + R), where utj = d^~^^u/dtdx^ , j = 0,r'. As a result, we obtain the following 
condition 

Ft + Fu^DiiSuxx + R) + Gx + Gu,uj+i = 0. (11) 

Let us decompose pT]) with respect to the highest derivatives uj. Thus, the coefficients of ^^'+2 
and Ur'+i give the equations Fu^, = 0, Gu^, + SFu^,__^ = that implies 

F = F, G = -SFu,,_,Ur'+G, 

where F and G are functions t, x, u, ui, . . . , Uj.i_i. Then, after selecting the terms containing 
ul,, we obtain that - S F^^, ^^u,., _^ = 0. It yields that F = F^u,j^i + where F^ and F° 
depend only on t, x, u, ui, . . . , Ur'-2- 

Consider the conserved vector with the density F = F — DxH and the flux G = G + D^H, 
where H = J F^dur'-2- This conserved vector is equivalent to the initial one, and 

F = F{t, X,U,Ui, . . . , Ur'-2), G = G{t, X,U,Ui, . . . , Ur'~l). 

Iterating the above procedure a necessary number of times, we result in an equivalent con- 
served vector depending only on t, x, u and Ux, i.e. we can assume at once that r' < 1. Then the 
coefficients of Uxxx and Uxx in ([H]) lead to the equations Fu^ = 0, Gu^ + SF^ = that implies 
F = F{t, X, u) and, moreover, G = —F^ J S dux + G^, where G^ = G^{t, x,u). □ 
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Note 5. A similar statement is true for an arbitrary (l+l)-dimensional evolution equation L 
of the even order r = 2f, f G N. For example [9], for any conservation law of L we can assume 
up to equivalence of conserved vectors that F and G depend only on t, x and derivatives of u 
with respect to x, and the maximal order of derivatives in F is not greater than f. 

Lemma [2] gives a stronger result for a more restricted class of equations. In the above proof 
we specially use the most direct method to demonstrate its effectiveness in quite general cases. 
This proof can be easily extended to other classes of (H-l)-dimensional evolution equations of 
even orders and some systems connected with evolution equations [26]. 

Theorem 2. A complete list of Gi -inequivalent equations ([2]) having nontrivial conservation 
laws is exhausted by the following ones 

1. h = l: ifu, -Au,- JB). 

2. A = l, f = -h{h~^),,: {e\h-^),,u, e\h-\,-{h'%u + jB)). 

3. B = eA + l, f = hy: {e*-'' f'hyU, -e^Ae'" f'^Uy + hu)). 

4. B = eA + l, f = hy + hy-^: {e^-" J'''yfu, -e^yAe'^ J'^Uy + yhu - JA)). 

5. B = eA: {yfe~^^u, -yAe'^-f^ + JA), {fe'^^^u, -Ae'^^^Uy). 

6. B = eA + l, f = -Z-\ h = Z-y^e.p[^-l^^^^dyy. 

((cr'=iy + o-^0)/e-^/^u, -{a'^^y + a''^){Ae'^ ■f'^Uy + hu) + a''^ JA). 

7. A = l, B = 0: {afu, —aux + axu)- 

Here y is implicitly determined by the formula x = J ^ ^^^^^'^^ dy ; e,aij = const, i,j = 0,1; 
{a^^,a^^) = {a'^^ (t) , a^^ (t)) , k = 1,2, is a fundamental solution system of the system of ODEs 
erf = a^ya^^; Z = ooiy^ + (aoo — — o-W, ct = cn{t,x) is an arbitrary solution of the linear 

equation fat + a^x = 0. Hereafter JA = J Adu, jB = J B du. In case [5] e G {0, 1} mod Gi . 

(Together with constraints on the parameter-functions A, B, f and h we also adduce conserved 
vectors of the basis elements of the corresponding space of conservation laws.) 

In Theorem [2] we have classified conservation laws with respect to the usual equivalence group 
Gi . The obtained result can be formulated in an implicit form only, and indeed case [6] is split 
into a number of inequivalent cases depending on values of aij. At the same time, using the 
extended equivalence group Gi , we can present the result of classification in a closed and simple 
form with a smaller number of inequivalent equations having nontrivial conservation laws. 

Theorem 3. A complete list of Gi -inequivalent equations ([2]) having nontrivial conservation 
laws is exhausted by the following ones 

1. h = l: ifu, -Aux-JB), 1. 

2. ^ = 1, S„/0, f = -h{h-%x- {e\h-%xu, e\h-^Ux-{h-^)xU + jB)), -e'h'\ 
3.5 = 1, f = hx: {e'fu, -e\Aux + hu)), e'. 

4. -6 = 1, f = hx + hx^^ : ( e^xfu, —e^{xAux + xhu — JA) ), e*x. 
5a. -6 = 0: {fu, —Aux), 1; (xfu, —xAux + f^), x. 

5b. -6 = 1, / = 1, h = l: {u, -Aux-u), I; {{x + t)u, -{x + t){Aux + u) + jA), x + t. 
5c. 5 = 1, / = e^, /i = e^: ( e^+^-u, -e*(^n^ + e^u) ), e*; 

(e^+*(x + t)n, -e\x + t){Aux + e''u) + e* jA), e*(x + t). 
5d.5 = l, f = xi'-^, h = x'': {e''^x''-^u, -e^'\Aux + x''u)), e''^; 

( e^'^+^^^xf'u, e(^+i)*(-x^Ua; - x^+^u + jA) ), e^^^+^^^x. 
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6. B = l, / = e-^/^a;-3, /i = e-/^/^x-\ ^uG {0,1}: 
{e^^xfu, —e'^^x{Aux + hu)+e'^^fA), e^^^x; 
{e^'\tx-l)fu, -ef'^itx -l){Aux + hu) +te^'^ JA), e^*(te-l). 

7. B = 1, f = \x - l|/^-3/2|x + l|-M-3/2^ h=\x- l\^'-^/^\x + Ij-A'-l/S . 

( e(2M+i)*(x - l)/u, -e(2^+i)*(x - l)(^u^ + /in) + e^^^+i)* /A ), e(2A'+i)* (x - 1); 
( e(2M-i)t (x + l)/n, -e(2/^^i)*(x + l)(^u^ + /in) + e^^^-i)* ), e(2/'-i)t (x + 1). 

g ^ _ J _ g/i arctan a; _j_ j^^— 3/2 ^ _ g/^ arctan x |^^2 _j_ -^^ — 1/2 . 

( e'^*(x COS t + sint)/n, —e'^*(x cost + sin t)(Ana; + /in) + e'^* cost JA), e'^* (x cos t + sin t) ; 

( e'^*(x sin t — cos t)fu, —e^^{x sin t — cos t)(^n2: + hu) + e'^* sint jA ), e^*(x sin t — cos t). 

9. A = l, B = 0: {afu, —aux + axu), a. 

Here fi = const, a = a{t,x) is an arbitrary solution of the linear equation fat + axx = 0. 
(Together with constraints on the parameter-functions A, B, f and h we also adduce conserved 
vectors and characteristics of the basis elements of the corresponding space of conservation laws.) 

Proof. In view of lemma [21 we can assume at once that F = F{t,x,u) and G = G{t, x,u,Ux). 
Let us substitute the expression for ut deduced from ([2]) into ([H) and decompose the obtained 
equation with respect to Uxx- The coefficient of Uxx gives the equation AFu+ fGu^ = 0, therefore 
G = —Af~^FuUx + G^{t,x,u). Taking into account the latter expression for G and splitting 
the rest of equation ([9|) with respect to the powers of Ux, we obtain the system of PDEs for the 
functions F and G^ of the form 



-bf^-a(^ 

f V / 

Solving first two equations of (fT2]) yields 



Fuu = 0, -BFu-A(^) +Gi = 0, Ft + Gl = 0. (12) 



F = F\t,x)u + F%t,x), G^ = (^) jA-t^F^ JB + G^{t,x). 

^ J J 

In further consideration the major role is played by a differential consequence of system (jl2p 
that can be written as 

Indeed, it is the unique classifying condition for this problem. In all classification cases we obtain 
the equation Fj) + = 0. Therefore, up to conserved vectors equivalence we can assume F^ = 
G^ = 0, and additionally F^ for conservation laws to be non-trivial. Equation (|13|) implies 
that there exist no non-trivial conservation laws in the general case. Let us classify the special 
values of the parameter- functions for which equation ([2|) possesses non-trivial conservation laws. 
There exist four different possibilities for values of A and B. 

1. dim{A,B, 1) = 3. It follows from ([13]) that F/ = {F^ / f)xx = {hF^ / f)x = and therefore 
F^ = C1///1, {\./h)xx = 0, i.e., obviously h G {l,x^'^} mod Moreover, /i = 1 ~ /i = x~^ 
mod G~ (the corresponding transformation is x = In |x| and / = x^/, the other variables and 
parameter-functions are not changed). As a result, we obtain case[TJ 

2. A G (1), B ^ (1). Then A = 1 mod G~ and {hF^/f)x = 0, F/ + {F^/f)xx = 0, i.e., 
F^ = a{t)f/h, where at/a = x. = const, x 7^ (otherwise we have case [T|) and so x = 1, 
/ = —h{h~^)xx mod G^ (case [2]). 
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3. ^ (1), B G Then B G (1) mod G~ and (F^f)^^ = 0, = B{hF^/f)^, i.e. 
= (a^(t)j; + a^{t))f and a|x/ + a^/ = i?(a-^(x/i)2: + a^hx). For S = we obtain case[5^ at 

once. Suppose B ^ 0. Then B = 1 mod G~ and the dimension m = dim{f,xf,hx, {xh)x) can 
have only the values 2 and 3. 

If ?Ti = 3 then there exist constants a^^, bfj_, /i, = 0, 1, and a function 6 = 6{x) such that 
{bo.hi) / (0,0), dim(/,x/,6') = 3 and = aoof + aoixf + boO, {xh)^ = a^f + auxf + bi6. 
Therefore, aj' = afj,ua^i bfj,a^ = 0, i.e. = C^e°"*, where C^,(T = const and cr 7^ (otherwise, 
this case is reduced to a subcase of d]), hence a = 1 mod Depending on values (either 
vanishing or non- vanishing) of Ci we obtain cases [3] and [5] correspondingly. 

If m = 2 then hx = aQof + oqix/, (x/i)a; = aio/ + auxf for some constants a^i,, fi,!^ = 0, 1. 
Therefore, = a^j^a^, / = —h/Z, hx/h = —{aoix+aoQ)/Z, where Z = aoix'^+(aoo—aii)x—aiQ, 
i.e., h = exp(— ^a^^ J Z~^(ix). As a results, we obtain two conservation laws with the 

conserved vectors 

{a'\t)x + a'°{t))fu, -{a'^{t)x + a'^{t)){Aux + hu) + o^^{t) /A), 

where (a*^, a**^), i = 1,2 form a fundamental set of solutions of the system = a^^a^. Separate 
consideration of possible inequivalent values of the constants a^^ leads to cases [5)3-[51i and [GHSl 

4. A,B (1). Therefore, A = 1, S = mod and + {F^ / f)xx = (caseE]). □ 

Analysis of the classification given in theorems [2] and [3] results in the following conclusions. 

There are no conservation laws for equations from class ([2]) in the general case. Imposing 
restrictions on values of arbitrary parameters leads only to subclasses of ([2]) with one-, two- 
and infinite-dimensional spaces of (local) conservation laws. An equation from class ([2]) has 
infinite number of linearly independent conservation laws iff it is linear. In this case the space 
of conservation laws is parameterized with an arbitrary solution of the corresponding backward 
equation (which is obtained from the initial one by the time reflection). 

There are four subclasses of equations ([2]) with one linearly independent conservation law 
(cases [iHll) ■ Case [U is distinguished from the other ones due to the single constraint h = const 
and, therefore, arbitrariness of nonlinearities A and B and the parameter-function /. The 
distinguishing features of case [2] are disappearance of nonlinearity with respect to A due to the 
constraint A = const, arbitrariness of B and a second-order differential constraint between / 
and h. All the latter cases are obtained with the constraint B G {A, 1) and first-order differential 
constraints of / and h since only cases [3] and H] with the constraint B G (^,1) admit extensions 
of space of conservation laws. 

Cases [Sb-Oi of Theorem [3] can be reduced to case [5^ by means of additional equivalence 
transformations which belong neither to G~ nor even to G~ and are pure form-preserving point 
transformations [20-22] for class ([2]): 

[5b ^r5b/-i : i=t, x = x + t, u = u; 
[5h -^~5kf-ex : t = e*, x = x + t, u = u; 

[5ll(^ + 1 / 0) ^3Jij=^M-i : t = (/i-M)"^(e(''+^)* - 1), i = xe\ u = u; 

[5li(;U + 1 = 0) ^l5^j=j,-2 : i=t, x = xe\ u = u. (14) 

6 Contractions of conservation laws 

Similarly to contractions of symmetries or equations [15] one can consider the notion of contrac- 
tions of conservation laws. 

Consider the class {C{e)} of systems C{e): L{x, U(^p),e) = of / differential equations for m un- 
known functions u = (n^, . . . , u"^) of n independent variables x = {xi, . . . , Xn), which are param- 
eterized with the parameter s. Here U(p) denotes the set of all the derivatives of u with respect to x 
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of order not greater than p, including u as the derivatives of the zero order. L = [L^, . . . , L') is a 
tuple of / fixed functions depending on x, and e. Let A{e) = {X^{x, e), . . . , X''{x, li(fc), e)) 
are characteristics of conservation laws of the systems C{e), i.e., A(e) € Ch(>C(e)). Suppose also 
that A(e) ^ A, C{e) ^ C, e ^ in C''{JP). 

Consider the action of Euler operator E on characteristic forms of the conservation law of 
the systems £{e): EX^'{e)L^'{e) = Ve. Since X^'{e)L'^{e) X^'L^',_e 0_in C''{JP) then 
= E{X''{£)L^'{e)) E{A^'L^') = 0, e ^ 0. Therefore E(A^L'^)_= or A G Ch{C) and A^Z^ = 
is a characteristic form of the conservation law of the system £.. 

By analogy with terminology accepted for Lie algebras, we will call such limits as contractions 
of characteristics and contractions of conservation laws. 

Example 1. Consider the equation 

X^^'^Ut = {AUx)x + X^Ux 

having two linearly independent conservation laws with the characteristics A^* = e'^* and A2 = 
^(fj.+i)t^ (CaseOd of Theorem[3]). Under the contraction x = 1 + t = /i ^ +00 it goes 
to the equation 

e^u^ = {Aux)x + e'^Ux 

which possesses linearly independent conservation laws with the characteristics Ai = e* and 
A2 = e\x + i) (CaseOc of Theorem [3D. 

Under the same limit process the characteristics A^* and Ag are transformed to characteristics 
of the target equation. Name A^ = e^"* = e* ^ e* = Ai and /i(A^-A^) = e*((e*/^- l)//+xe*/'') 
A2 = e\i + i). 



Example 2. 



X — 1 



x+1 



l™2 



X — 1 



x+1 



-11-1/2, 



Under contraction x = 2ij,x/ij,', t = /i't/(2;u), fi +00 it maps to 

--Sg-Aj'/a;^, = {AUx)x + EX'^ e~^' Ux 

The characteristics A^* = e^'^^~^^^^{x — 1) and Ag = e^'^^~^^^{x + 1) of the conservation laws of the 
origin equation give rise the characteristics for the target equation: 



= e(2/^+i)*(x - 1) = e^'''—-- ( X - ^ ) - e^'*x = Ai 



and 



i(Ai - A2) = e^'* ( x ^'''^' ^ - i(e^'*/^ + 1)) e^'~"^^^^ - e^'\ii - 1) = A2. 
2 \^ ^'/^i 2 J 

The problem of finding all contractions of conservation laws of class ([2]) is closely connected 
to the problem of construction of contractions of equations ([2]), and therefore, still remains open. 
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7 Generating sets of conservation laws 

of nonlinear diffusion— convection equations 



We emphasize that the conservation laws adduced in Theorem [3] are G~-inequivalent and hnearly 
independent. At the same time, acting by the Lie symmetry transformations of equations having 
multidimensional spaces of conservation laws or by the equivalence transformations of the whole 
class of equations, one can find generating sets of conservation laws with respect to the symmetry 
groups. 

Note 6. Since the cases [Sb^Sli of Theorem [3] can be reduced to case [5^ by means of point 
transformations (jl4p . below we exclude them from the consideration and assume them to be 
equivalent to caseE^. 

Note 7. Below we investigate potential systems and potential conservation laws of nonlinear 
equations only. The detailed analysis of linear case can be found in [27] . 

Example 3. It is obvious that in general, for case 5a there exists no symmetry transformation 
changing the conservation laws of the fixed equation. However, for subcase / = 1 of case 5a 
the space translation x ^ x + 1 maps the second basis conserved vector to ( (x + l)u, — (x + 
l)Aux + JA)- Subtracting it from the second basis vector, we obtain {—u,—Aux), that is 
the first basis vector multiplied by —1. Therefore, the generating with respect to symmetry 
group set of conservation laws of the subclass / = 1, S = of equations ([2]) consists of the 
Dt{xu) + DxijA - xAux) = 0. 

Investigation of generating sets of conservation laws can be performed both in terms of 
conserved vectors (as in above example) or in terms of characteristics. Below for short of 
presentation we use the characteristics terminology. 

Example 4. Consider case [6] of equations possessing two-dimensional space of conservation 
laws with the basis conserved vectors having characteristics e'^'x and e'^*(tx — 1). Symmetry 
transformation of time translation t — > t + 1 maps the second characteristic (modulo multiplying 
by a constant) to e^^{tx — x — 1). Subtracting the result from the characteristic e'^*(tx — 1) 
we get exactly the first basis characteristic e'^^x. Therefore, the generating set of conservation 
laws of the subclass 6=1,/ = e~^^^x~^, g = 1, h = e~^^^x~^ of equations ([2]) consists of the 
conservation law with characteristic e^^{tx — 1). 

Similarly one can prove that conservation laws in case 8 are generated by one with the 
characteristic e^^{xcost + sint) and in case 7 by the conservation law with the characteristic 
e^'^*(xcoshi — sinht). Note that the basis conservation laws in this case have the characteris- 
tics e^'^^~^^^{x + 1) and e^'^^^^^*(x — l)fu. None of these two characteristics generates the whole 
space of conservation laws. 

Let us consider now the problem of construction of generating set of conservation laws with 
respect to the extended equivalence group . 

Example 5. Consider subcase [5^ {B = 0) having two-dimensional space of conservation laws 
spanned by ones with characteristics 1 and x. Subtracting the second characteristic from the 
result of application of equivalence transformation x = x -|- 1 to it, we get the first basis con- 
servation law for equation fut = {Aux)x with /(x) = /(x — 1). Therefore, we can assume 
that conservation laws of subclass B = are generated by Dt{xfu) — Dx{xAux — JA) = and 
group G~ of equivalence transformations. 

Example 6. Let us investigate in more detail case 7 of Theorem [3l It is shown above that for 
the fixed value of fx the set of local conservation laws is generated (with respect to the symmetry 
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group of the equation) by the conservation law with characteristic e^^*(xcosht— sinht). Alterna- 
tively, we can consider an action of equivalence transformations to the conservation laws. Thus, 
discrete equivalence transformation i = —t, x = —x, fj, = —fi applied to the first basis character- 
istic e^'^^^^^^ {x — 1) generates the second characteristic e^^^""*^-'* of equation with / = — /(— x, —fi) 
h = —h{—x, — /i), A = —A and jl = —fi. 

Similarly it can be shown that there exists no other nontrivial action of transformations 
from G~ to conservation laws of equations ([2]). Therefore, the following theorem holds. 

Theorem 4. The generating with respect to set of nonlinear equations ^ and corresponding 
conservation laws consists of 



1. 


h-. 


= 1: 


{fu, -Au^-JB); 


2. 


A 


= 1, 


Bu^O, f = -h{h-^),,: (e*(/i-i),,n, e*(/i-iu. 


3. 


B 


= 1, 


f = hj;: [e^fu, -e\Aux + hu)); 


4. 


B 


= 1, 


f = hx + hx^^: {e^xfu, —e^{xAux + xhu — jA)) 


5. 


B 


= 0: 


(xfu, -xAux + JA); 


6. 


B 


= 1, 


/ = e-z^/^x-s, h = e-^'/''x-\ fx e {0,1}: 



{e''\tx - l)fu, -e^'\tx - l){Aux + hu) + te^"^ jA); 
7. B = 1, f = \x - l\^'-^/'^\x + l|-/^-3/2^ h=\x- l\^'-^/'^\x + l|-^-V2 ; 
(g(2M+i)t(^ - -e(2M+i)*(2; - l){Aux + hu) + e(2M+i)t jA); 

g — Y f = a-rctan a; |^^2 _j_ j^^— 3/2 ^ _ g/^ arctan x ^^^2 _j_ X)~^/2 • 

(e'^*(xcost -I- sint)/u, —e'^*{x cost + smt){Aux + hu) + e^* cost jA). 

8 Potential systems and potential conservation laws 

If the local conservation laws of a system C of differential equations are known, we may apply 
Lemma [1] to the constructed conservation laws on the set of solutions oi C = . In such way 
we introduce potentials as additional dependent variables. Then we attach the equations con- 
necting the potentials with components of corresponding conserved vectors to . (If n > 2 the 
attached equations of such kind form an underdetermined system with respect to the potentials. 
Therefore, we can also attach gauge conditions on the potentials to .) 

We have to use linearly independent conservation laws since otherwise the introduced poten- 
tials will be dependent in the following sense: there exists a linear combination of the potential 
tuples, which is, for some r' G N, a tuple of functions of x and U(r') only. 

Then we exclude the unnecessary equations (i.e., the equations that are dependent on equa- 
tions of and attached equations simultaneously) from the extended (potential) system 
which will be called a potential system of the first level. Any conservation law of is a one 
of C^. We iterate the above procedure for to find its conservation laws which are linearly 
independent with ones from the previous iteration and will be called potential conservation laws 
of the first level. 

We make iterations until it is possible (i.e., the iteration procedure has to be stopped if all 
the conservation laws of a potential system C^~^^ of the {k + l)-th level are linearly dependent 
with the ones of C^) or construct infinite chains of conservation laws by means of induction. A 
such way may yield purely potential conservation laws of the initial system C, which are linearly 
independent with local conservation laws and depend explicitly on potential variables. 

Any conservation law from the previous step of iteration procedure will be a conservation 
law for the next step and vice versa, conservation laws which are obtained on the next step 



15 



and depend only on variables of the previous step are linearly dependent with conservation laws 
from the previous step. It is also obvious that the conservation laws used for construction of a 
potential system of the next level are trivial on the manifold of this system. 

Since gauge conditions on potentials can be chosen in many different ways, exhaustive real- 
ization of the above iteration procedure is improbable in case n > 2. 

The case of two independent variables is singular with respect to possible (constant) in- 
determinacy after introduction of potentials and high effectiveness of application of potential 
symmetries. That is why we consider some notions connected with conservation laws in this 
case separately. We denote independent variables as t (the time variable) and x (the space one) . 
Any local conservation law has the form 

DtF{t, X, U(.r) ) + DxG{t,x, U(r) ) = 0. 

It allows us to introduce the new dependent (potential) variable v by means of the equations 

Vx = F, vt = -G. (15) 

In the case of single equation C, equations of form (jlSp combine into the complete potential 
system since C \s a. differential consequence of (I15p . As a rule, systems of such kind admit a 
number of nontrivial symmetries and so they are of a great interest. 

Lemma [3] and equation (jl5p imply the following statement. 

Proposition 4. [26] Any point transformation connecting two systems C and L of PDEs with 
two independent variables generates a one-to-one mapping between the sets of potential systems, 
which correspond to L and C. Generation is made via trivial prolongation on the space of 
introduced potential variables, i.e. we can assume that the potentials are not transformed. 

Corollary 2. The Lie symmetry group of a system C of differential equations generates an 
equivalence group on the set of potential systems corresponding to C 

Corollary 3. Let C\s be the set of all potential systems constructed for systems from the 
class C\s with their conservation laws. Action of transformations from G""{L,S) together with 
the equivalence relation of potentials naturally generates an equivalence relation on C\s. 

Note 8. Proposition H] and its corollaries imply that the equivalence group for a class of systems 
or the symmetry group for single system can be prolonged to potential variables for any step of 
the direct iteration procedure. It is natural the prolonged equivalence groups are used to classify 
possible conservation laws and potential systems in each iteration. Additional equivalences which 
exist in some subclasses of the class or arise after introducing potential variables can be used 
for further analysis of connections between conservation laws. 

9 Potential systems of diffusion— convection equations 

At first we consider potential systems obtained with introducing potentials associated to a single 
conservation law. According to the terminology proposed in [26] we use the attribute "simplest" 
to single out such systems from the set of all potential systems associated to equations . 

The spaces of local conservation laws in cases 1-4 of Theorem [3] are one-dimensional and give 
rise to the following nonlinear simplest potential systems: 

1. h = l: Vx = fu, vt = Aux + J B 

2. A = l, B^^O, f = -h{h-^)xx: Vx = e\h~^)xxu, vt = e\-h-^Ux + {h~^)xU - jB). 

3. B = 1, f = hx- Vx = e^hxU, vt = e^{Aux + hu). 
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4. B = 1, f = hx + hx ^: = e^xfu, = e^{xAux + xhu — JA). 

If the dimension of the space of local conservation laws is greater than one then the situa- 
tion is more complicated. Previously, for construction of simplest potential systems only basis 
conservation laws were used (e.g., like systems 1 and 5 for case 5a of Theorem [3|). It was just 
luck that all possible inequivalent potential systems are exhausted for the considered cases. It 
happened because arbitrary linear combination of basis conservation laws could be reduced to a 
subset of the basis elements by means of a symmetry transformation of the considered system. 
Here we will show that the basis conservation laws may be equivalent with respect to group of 
symmetry transformations, or vice versa, the number of C-independent linear combinations 
may be grater than the dimension of the space of conservation laws. Class ([2]) is very rich in the 
sense that it provides us with examples of all these three cases. 

Case 5a of Theorem [3] is the classical. The most general form of the simplest system is 

Vx = (CIX + C2)fu, Vt = (CIX + C2)AU;j: - J A, 

where ci and C2 are arbitrary constants. Here and below (ci,C2) 7^ (0,0). Moreover, without 
lost of generality, we can assume that the coefficients of linear combinations are determined up 
to nonzero factor (in most of cases cf + = 1). 

If ci = the obtained system coincides with system 1|_b=o. If ci 7^ it can be mapped to 
the system 

5. B = 0: Vx = xfu, Vt = xAux — J ^ 

by means of equivalence transformation of translation in x. We emphasize that in case /x 7^ 
this transformation changes the value of arbitrary element /. Therefore, we can reconstruct 
all potential symmetries of the subclass B = of class ^ applying translations in the space 
variable to the arbitrary element / and to the potential symmetries obtained from systems 1 
and 5. 

Potential system corresponding to case 6 has the form 

Vx = e'^*(cix + C2tx — C2)fu, Vt = e^*[cix + C2tx — C2){Aux + hu) — e^*ci jA. 

In contrast to the previous case, we can simplify the form of the potential system using the 
symmetry transformations (translation of time) of the fixed equation of the given form. More 
precisely, using the time translation and scaling of potential v one can reduce the potential 
system to one of the following forms: 

6.1. B = l, f = e-^/^x-^ h = e-^l^'x-^, G {0, 1}: 

Vx = e^^^xfu, Vt = ef^^x{Aux + hu) — e^* JA, 
if C2 = or 

6.2. B = 1, f = e~^/^x-^ h = e-''/^x-^ /i G {0, 1}: 

Vx = e^'^tx - l)fu, Vt = e^*(te - l){Aux + hu) - te^"^ /A 
if C2 ^ 0. 

The most interesting from this point of view is the case 7 of Theorem [3l It is obvious that 
the most general form of the associated simplest potential system is 

Vx = e^^*(x(ci cosht + C2 sinht) — ci sinht — C2 cosht)/'u, 
Vt = e^^*(x(ci cosht + C2 sinhf) — ci sinhi — C2 cosh. t){Aux + hu) 
- e^^* (ci cosh t + C2 sinh t) jA). 



17 



Without lost of generality we can assume ci > 0. Cases ci > |c2| > are equivalent (under the 
action of symmetry transformations of translation in time) to the system 

7.1. B = l, f =\X- l|^-3/2|x + l|-M-3/2^ h=\x- + 

Vx = e^^^* (x cosh t — sinht)/u, vt = e^^^^{xcosh.t — smh.t){Aux + hu) — e^'^^ cosht JA. 
Similarly cases |c2[ > ci > fall to the system 

7.2. B = l, f=\x- 1|'^"3/2|^ j_ l|-M-3/2^ h=\x- ll'^-l/^lx + 11-^-^/2. 

Vx = e^^^* (x sinh t — cosht)/u, vt = e^^^^{xsinht — cosht)(j4ua; + hu) — e^'^^sinht fA. 

At last, we should consider the case of [ci[ = \c2\- Using the scaling transformation of the 
potential v we can reduce them either to 

7.3. B = 1, f = \x - l|'^-3/2[^ ^ i|-M-3/2^ h=\x- l\^'-^/^\x + Ij-/^-^/^. 

= e(2M+i)*(x - l)/u, Vt = e(2/^+i)*(x - l){Aux + hu) - e^^M+i)* jA, 
if ci = C2 or to 

B = l, f =\x- l|^-3/2|^ ^ l|-/^-3/2^ /l = |x - ^ l|-M-l/2. 

= e(2M-i)*(x + l)/u, Vt = e(2^-i)*(a; + l)(^n^ + hu) - e^^^-i)* 

in the case ci = — C2. Since we consider simultaneously the whole class of equations ([2]), we 
can apply additionally discrete equivalence transformation of alternating of sign in the set 
{t, X, u, V, f, h, A, fi) and reduce the latter system to 7.3. We emphasize once more that if one 
considers a separate fixed equation from case 7 of Theorem [31 he obtains 4 independent poten- 
tial systems. Considering simultaneously the class of equations we have three G^^-independent 
systems. Nevertheless, in this case the number of independent potential systems is greater than 
the number of basis conservation laws. 

It is not difficult to show that under the action of symmetry transformation of translation in 
time and scaling of potential variable there exist only one locally inequivalent simplest potential 
system in case 8 of Theorem [3l 

8. B = l, f = e'^^'-'=t^'i^'(x2 + l)-3/2^ h = e/^^=-'=t^'i^(2;2 + l)-i/2. 

Vx = e^*(xcost + smt)fu, vt = e^^{xcost + sm.t){Aux + hu) — e^* cost jA. 

Thus, the number of independent potential systems is smaller then the number of basis conser- 
vation laws. 

Theorem 5. The complete set of Gi -independent simplest potential systems of equations from 
class ([2]) is exhausted by ones with bold numbers 1-5, 6.1, 6.2, 7.1-7.3, 8. 

As one can see, in cases 5a-8 of Theorem [3] the spaces of local conservation laws are two- 
dimensional. It allows us to introduce extended potential systems [26] by means of introducing 
two potentials for each case simultaneously: 

5'. B = 0: Vx = fu, Vt = Aux, Wx = xfu, wt = xAux — JA. 

6'. B = l, f = e-f/''x-^, h = e-'^/^'x-i, fi G {0, 1}: 

Vx = e^^xfu, Vt = ef^*x{Aux + hu) — e^* jA 

Wx = e''*(tx - l)fu, Wt = ef'^itx - l){Aux + hu) - te'"* JA. 
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7'. B = 1, f = \x - l|'^-3/2|a; + l|-M-3/2^ h=\x- l\f'-^/'^\x + lI-^'-V2. 

t,^ = e(2A'+i)*(x - 1)/m, r;^ = e^^^'+^^\x - l)(An^ + hu) - e^^M+i)* 

u;^ = e(2/^-^)*(x + u;j = e^^f'-^^^x + l)(An^ + hu) - e^^^-^)* p. 

8'. B=l, f = eA'arctanx^^2 _^ l)-3/2^ /j ^ g/. arctan x ^^2 _^ l)-l/2. 

= e^^{xcost + smt)fu, vt = e^*(xcost + sint)(Ana; + hu) — e^* cost /A, 
Wx = e^^{xsmt — cost)fu, wt = e^^{xs'mt — cost){Aux + hu) — e^*sint JA. 

Potential symmetries of equations ^ arising from the above independent potential symme- 
tries will be investigated in the last part [16] of this series. 

10 Potential conservation laws 

of diffusion— convection equations 

All potential systems of equations ([2]) are constructed with usage of local conservation laws of 
equations ([2]) which were classified with respect to . Each of these subclasses of conservation 
laws is equivalent with respect to a subgroup of Gi . In each case this subgroup has a very 
simple structure and can be singled out from ^ imposing an additional condition = 5s = 0. 
It becomes the usual equivalence group of class ([2]) and is trivially prolonged to the corresponding 
potentials. Henceforth we will denote such prolongation as G~ . 

Let us investigate local conservation laws of potential systems 1-8', which have the form 

DtF{t,X,U(^j.),V(_r){,W(^r))) + D^G{t,X,U(^^),V(^r){,W(^rj)) = 0. (16) 

These laws can be considered as nonlocal (potential) conservation laws of equations from 
class 

Lemma 3. Any conservation law of form ()16p for each of systems 1-8' from Table 1 has the 
zero order, i.e., it is equivalent to a conservation law with a conserved density F and a conserved 
flux G that are independent on the (non-zero order) derivatives of u and potentials v (and w). 

Proof. Consider any from the systems 1-8'. Taking it and its differential consequences into 
account, we can exclude dependence of F and G on the all (non-zero order) derivatives of v (and 
w) and the derivatives of u containing differentiation with respect to t. The remain part of the 
proof is completely similar to the one of Lemma [2j □ 

Systems 1|_b=o, 3-7.3 and 5'-8' are quit similar: for all of them B = const. Analyzing 
conservation laws of them we prove the following statement. 

Lemma 4. Nonlinear equations ^ with B = or B = 1 have nontrivial potential conservation 
laws only if A = n"^ mod Gpj, . 

Proof. Each of the simplest potential systems l|5=o, 3-7.3 has the form 

Vx = Xfu, Vt = X{Aux + hu) - Xx J A, (17) 

where A = X{t,x), is a characteristic of a local conservation law of equation ([2]) with B = const. 
According to Lemma [3] we look for its conservation laws in form 

DtF{t, x, u, v) + DxG{t, X, u, v) = 0. 
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First we expand this expression on the solution manifold of (jl7p and decompose it with respect 
to unconstrained derivatives of u. Coefficients of ut and Ux give 

Fu = 0, G = -XF^JA + Git,x,v). 

Substituting this to the rest of the conservation law we obtain a classifying equation for the 
components of the conserved vector: 

Ft + F,{Xhu -X^JA)- {XF,)x JA + Gx- X^fF,,u j A + XfG^u = 0. (18) 

Uuf A^ {l,uj A) then 

Fvv = 0, -Ax-K - {XF^)x = 0, Fyh + fG^ = 0, Ft + G^ = 0. 

Therefore, 

F = F\t, x)v + x), G = -jF^v + x). 

Up to the usual equivalence of conserved vectors we can assume that F^{t,x) = G^(t,x) = 0. 
Consider an equivalent conserved vector 

F = F- Dx{^{t, x)v), G = G + Dt{^(t, x)v), 

where <^a;(t,x) = F^. It is obvious that this conservation law is local for the equation and 
moreover, it is trivial on the solution manifold of the potential system (]17p . 

I{ u J A G {1, u, f A) , then u J A = co + c\u + C2 / A. This is an algebraic equation for u j A. 
Up to translations of u we can assume that C2 = 0. Therefore u / A = cq + c\u ox A a {l,u^^} 
mod Gpj.. Since we consider only nonlinear equations, A = u~'^ is the only case where potential 
conservation laws can exist. 

Consider now the general potential systems 5'-8' on two potentials. Each of them has the 
form 

vi = X'fu, vi = X\Aux + hu)-Xl.jA, i = 1, 2, (19) 

where A* = A*(t, x) are linearly independent characteristics of local conservation laws. Similarly 
to the case of simplest potential system we can prove that an arbitrary conservation law of 
system (jl9p has the form 

DtF{t,x,v^,v'^) + D,^G{t,x,u,v^,v'^) = 0, 

where the flux G = —X'F^i J A + G{t, x, f^, u^), and the classifying equation is the following 

Ft + G, + uX'{hF,. + fG,^) - J A{XIF,. + {X'F,.),) - X'X^fF,^,,u J A = 0. (20) 

Now, ii u J A ^ {l,u, J A) (in other words, A ^ 1, mod Gpj.) we can split (llSp with respect 
to the linearly independent functions of u and obtain 

Ft + Gx = 0, X\hF,. + /G.O = 0, XiF,. + (A*F,0. = 0, X'X^fF,.,, = 0. 

The latter system implies, in particular, that 

F = F^{t,x,uj)e + F^{t,x,uj), G = -jF^e + G°{t,x,u;), 

2F'X^Xi + X^X^F^ixy - XW) = -F^iX^Xl - X^Xl)e - F^X^Xl - F^^{X^Xl - A^A^), 
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where uj = X^v'^ — X^v^, 9 = X^v^ + A^w^. Decomposing {X^v^ — X^v^) with respect to uj, 6 and 
collecting the coefficients of 9 in the last equation we obtain 2F^(A^A^ — A^A;j.) = 0. One can 
easily check that for each of the systems 5'-8' A^A^ — A^A;|. / 0, therefore, -F^^ = which implies 
immediately F^^ = 0. Thus, we get that the conserved density and, consequently, the conserved 
flux are linear with respect to the potentials. Therefore, in view of [6,25], the conservation law 
is local for equation ([2]). 

The lemma is completely proved. □ 

In view of Lemma HI in order to construct potential conservation laws of equations ([2]), it 
is enough to investigate conservation laws of the potential systems 1, 2 and cases A = u~'^ of 
systems 3-8'. The chain of lemmas below describes completely the set of potential conservation 
laws of equations ([2]) . 

Lemma 5. The list of C^^-inequivalent linearly independent conservation laws of system 1 is 
exhausted by the following ones 

1. VA, fB = ujA, f = l: Dt{e-) + D^{e- J A) = 0, 

2. A = u-\ B = 0, / = 1 : Dt{a) + D^{a^u-^) = 0, 

3. A = 1, B = 2u, f = 1 : Dtiae") + D^{a^e'' - aue") = 0, 

where a = a{t, x) and a = cr(t, v) are arbitrary solutions of the backward linear heat equations 
o-t + otxx = and at + (Tw = 0, correspondingly. 

Proof. In view of Lemma [3] we look for conservation laws of system 1 in form 

DtF{t, X, u, v) + DxG{t, X, u, v) = 0. 

Expanding the total derivatives in the expression for the conservation law and decomposing it 
with respect to unconstrained derivatives of u we get 

= 0, G = -F,fA + Git,x,v), 
Ft + F^jB-F^xjA + Gx- fF^^u JA + XfG^u = 0. 

To integrate the latter equation one should consider separately three cases: 

• l,u, J A,u J A, J B are linearly independent; 

• l,u, f A,u f A are linearly independent and l,u, J A,u J A, J B are linearly dependent; 

• uJ AG {l,u,J A). 

All conservation laws in the first case are trivial. In the second case we have J B = cq + ciu + 
{C2U + C3) J A. Substituting this into the classifying condition and solving it up to G^^. we obtain 
case 1 of the list of the lemma statement. At last, ii u J A € {l,u, J A), then similarly to 
LemmaUl A E {1,^"^} mod G^^. that leads after some obvious calculations to cases 2 and 3 of 
the lemma statement. □ 

Lemma 6. System 2 has no nontrivial conservation laws. 

Proof. The proof is completely similar to the proof of Lemma O □ 

Lemma 7. Systems 3, 4, 6.1-8 have no nontrivial conservation laws. There exists exactly one 
nontrivial conservation law for system 5: 

A = B = 0, / = : Dt{x-'^cr) + Dxix'^cr^u-^) = 0, 

where a = a(t, v) runs the solution set of the backward heat equation at + a^y = 0. 
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Proof. Each of the systems 3-8 has the form (|17p . The conserved density and flux of the 
conservation law has the form 

F = F{t,x,v), G = -XF^fA + G{t,x,v). 

In view of Lemma [Hit is enough to consider case A = u^^ only. Then (jlSp implies 

X^F^ + (AF„)^. = 0, XhF^ + XfG^ = 0, Ft + G^ + X^fF^^ = 0. 

Therefore, up to the equivalence relation of the conserved vectors F = X^^a, G = — jA~^o", 
where a = a{t, v) satisfies the classifying equation 

(A-V)t - a + fa,, = 0. 

Integrating this equation for the given values of characteristics A implies the lemma statement. 

□ 

Lemma 8. The list of G'^^-inequivalent linearly independent conservation laws of system 5' is 
exhausted by the following ones 

1. A = B = 0, f = l: Dt{a^) + D,{al^u-^) = 0, 

2. A = B = 0, f = x-^: Dt{x-^a^) + Z)^(x- V^^^-i) = 0, , 

0"' = a^{t,v^), i = 1,2 run the solution sets of the backward heat equations a] + 0"*,^^ = 
correspondingly. Systems 6' -8' have no nontrivial conservation laws. 

Note 9. These potential conservation laws of equations ([2]) indeed coincide with conservation 
laws of the simplest potential systems 1 and 5. 

Proof. Substituting A = into the classifying condition (fTSl) and splitting it with respect to 
the linearly independent functions of u we obtain the following system of determining equations 
for the coefficients of the conserved vectors of the conservation laws of systems (I19p : 

X\hF,. + fG,.) = 0, A^F,. + (A^F.Ox = 0, Ft + G^ + XX^fF,.,^ = 0. (21) 

From the first equation of ()2ip we obtain that G = - jF + H{t, x,u)). Here, as in Lemma HJ 
uj = X^v^-X^v\e = X^v^ + X^v^. 
It follows from (UHl) that 



X\l - X^vl = = u;, - (A>2 _ xlv^), 

x\l - x\] = {x'xl - x'xDu-' =ut- {\W - A^'). 

Decomposing X^v"^ — A^t;^ and X\v'^ — X'fv^ with respect to {uj,0) from the latter equations we 
get 

\k\k \1\2 \2\1 \k\k \1\2 \2\1 \1\2 \2\1 

^x^ . ^ '^x ~ ^ '^x n . \^ , ^ \ - ^ \ a , ^ '^x ~ ^ '^x 



Ux = -T^t^ T^T" = T-T-'^ + 



XJXi X^X3 X^X^ X^X^ u 

These equalities allow us to find the simplest form of conserved vectors corresponding to the 
conservation law. For now we have F = F(t,x,u:,9), G = X^ X^ Fgu~^ — jF + H{t,x,uj). An 
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equivalent conserved vector can be written as 

F = F + D^^{t, x,uj) = F + ^^ + ^^io^, 

G = G- Dt^{t,x,uj) = X^X^Fgu'^ -jF + H-^t- ^^ivt 

X 1 , 1 _i h ~ T „ / h A^A^ — A^AI A^A? — A^AJ 
X^X^Fqu ^ - -F + <^^e ( ^-r-p — -+ ' * 



/ V / A-^AJ AJAJ 

h ( X^X^ \ X^X^ 

/V AJA-' J XI XI 

For all of the characteristics corresponding to systems 5'-8', coefficients of in the expression 
for G is vanishing. Therefore, choosing ^ in such way that 

H + -A^x + <^..^-^') - - ^^^^ = 0, 



/V XI XI J " A^AJ 

we obtain that G = X^ X^Fgu^^ — jF, i.e., H = 0. That is why, without lost of generality, we 
assume that 

F = F{t,x,v\v'^), G = X'X'Feu-^ ~ ^ X'F^.u'^ - jF, 
where the conserved density F satisfies the following system: 

A^F,. + (A%), = 0, X'X^F,.,, + ^-j(^jf'^ =0. 
Taking into account its nontrivial differential consequence 

we can write the compatible system of determining equations for the conserved density F: 
^-U^f) =K{t,v\v^), X'X^F,.,,+K = 0, 



f f\f 

A;F,. + (A%), = 0. (22) 

System 5' is distinguished from the set of the systems (jl9p since we can assume h = 0. Let 
us investigate it separately. Substituting h = and the characteristics A^ = 1, A^ = x into ()22p 
and taking its differential consequence, one can obtain 

f,K,i + (xf, + 2f)K,2 = 0. (23) 



There exist three different cases of integration of equation (j23p : / is arbitrary function such that 
fx and xfx + 2/ are linearly independent, fx = and xfx + 2/ = cfx, where c is a constant. 

If fx and xfx + 2/ are linearly independent then K^i = K^2 = 0. Then, it follows from (|22p 
that Kt = 0. Therefore, K isa constant, and without loss of generality we can assume K G {0, 1}. 
If ii: = 1, then ([221) implies 



F = ft + M{x,v\v'^), M = -^9^ + M\x,lo)9 + M^{x,uj), 
iX'M^^)x + XiM^^=0. 
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Collecting coefficients of 9 in the last equation we get —x6 = 0, that is contradiction. Therefore 
K = Q. Then F^i = G^i = 0, and the conservation law is local. 

li fx = then / = 1 mod K^2 = 0. Then, from (j22p we obtain K^i^i + = and F = 
a^{t,v^) + M{x,v^ ^v"^), where K = a} = — cr^i^i. Substituting it into the second equation of ([22]) 
leads to the equation for M of form X^X^M^i^j = 0. Therefore, M = M^{x, uj)9 + M°(x, uo), and 
M^, satisfy the equation 

M° - xujMI + mI{x^ + 1) + 3xM^ = 0. 

Considering equivalent conservation law with the conserves vector F = F — Dx{N9), G = 
G + Df{N9), where N^{x,uj) = IVF one can show that iV = 0. Thus we get the conservation 
law from case 1 of the lemma statement. 

If xfx + 2/ = cfxi then up to translations of x we can assume that c = and / = x~'^. The 
local transformation i = x = —1/x, u = xu, = v^, v'^ = reduces this problem to the 
previous case / = 1. 

Case /i 7^ (systems 6'-8') can be investigated in a similar manner. At first, we find the 
differential consequence of system (I22p : 

(A7x. + 2A1./)K,. + F,. + £ (y)^ + 2AL - 2ALy) 

+F,.x (2y{j)^-xij + xi^=o. 

Considering the given sets of characteristics for each of the systems 6'-8' we obtain 

iyfx + 2Xif)K,.=0, 

where (A*/x + 2A^/) are linearly independent. Therefore, Kyi = Kf = 0. Absolutely similarly 
to case h = we get that the conservation law is local. 

The lemma is proved. □ 

As one can see, equation ([2]) admits nontrivial potential conservation laws of the first level iff 
it is linearizable or satisfies the conditions f = h = l, jB = ujA. It is shown in [26] that such 
equations do not have potential conservation laws of the second level. Therefore, the following 
lemma is true. 

Lemma 9. Any equation of form ([2]) has no nontrivial second level potential conservation laws. 

This statement completes investigation of potential conservation laws of the variable coeffi- 
cient diffusion-convection equations ([2]). Summarizing the above results, we can formulate the 
following theorem. 

Theorem 6. A complete set of G^^-inequivalent potential conservation laws of nonlinear equa- 
tions ([2]) consists of the following ones: 

1. ^A, jB = uJ A, f = 1: Dtie^ + Dx{e- f A) = 0, 

2. A = u'"^, B = 0, f = 1: Dt{cr) + DxiayU'^) = 0, 

3. A=l, B = 2u, f = l: Dt{ae") + Dxiaxc" - aue'') = 0, 

where the potential variable v{t,x) satisfies the potential system Vx = u, Vt = Aux + J B, 

4. A = u-^, B = 0, f = x~^ : + = 0, 

where Vx = x~^u, vt = xu~'^Ux + u~^. Here a = a{t,x) and a = cr(t,v) are arbitrary solutions 
of the backward heat equations at + Oxx = and at + (7yy = 0, correspondingly, 



24 



11 Conclusion 

In this paper we generalized the notion of equivalence of conservation laws with respect to groups 
of transformations and adduced description of generating sets of conservation laws. This frame- 
work is applied to investigation of conservation laws of variable coefficient diffusion-convection 
equations. We classified the local conservation laws of equations ([1]) both up to the usual group 
of equivalence transformations and up to the extended group containing nonlocal (with respect 
to the arbitrary elements) transformations. Usage of equivalence of conservation laws with re- 
spect to the extended equivalence group and correct choice of gauge coefficients of equations 
allow us to obtain clear formulation of the final results. 

The notions of contractions conservation laws and ones of characteristics of conservation laws 
are introduced. We presented some examples of such contractions in class ([T|). The problem of 
finding all possible contractions of equations and conservation laws in class ([T]) remains open. 

The proposed view on equivalence of conservation laws give rise to a generalization of the 
theory of potential symmetries. Namely we showed how to find all possible inequivalent potential 
systems associated to the given system of differential equations. We classified such potential 
systems and potential conservation laws for equations ([T|). 

Using this classification one can find all potential symmetries of class ([I]). This is a subject 
of the sequel part [16] of this series of papers. 
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